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Предлагается использовать метод реконструктивной томографии для расчета распределения зарядов и емкостей плоских про- 
водников неканонической формы. Приводятся примеры восстановления распределения зарядов и ёмкостей на модельных за- 
дачах. Сравниваются результаты известных расчетов ёмкостей с томографическим подходом. 


В реальных электротехнических устройствах 
происходят сложные процессы, связанные с нако- 
плением, перераспределением энергии, протека- 
нием токов смещения, проводимости. Поэтому 
при проектировании линий электропередач, кабе- 
лей, различных электрических и радиотехнических 
устройств возникает необходимость расчета элек- 
трической емкости. 

Известно, что емкость зависит только от геоме- 
трических параметров проводников и диэлектриче- 
ской проницаемости окружающей среды. В научно- 
технической литературе достаточно подробно рас- 
смотрены основные методы расчета емкости для 
проводников канонической формы, существует 
справочная литература по расчету емкостей [1], со- 
держащая примеры расчета, готовые формулы, та- 
блицы, численные результаты для проводников ка- 
нонической формы. Для проводников более слож- 
ной формы в работах [1-3] применяют специаль- 
ные методы, дающие приближенные результаты. 

Как альтернативу существующим методам рас- 
чета емкостей плоских проводников авторы пред- 
лагают метод компьютерной томографии, претен- 
дующий на высокую точность расчетов даже для 
проводников сколь угодно сложной конфигурации. 

Идея компьютерной томографии состоит в вос- 
становлении (реконструкции) функции двух пере- 
менных о{х,у) (распределения заряда) по их инте- 
гральным характеристикам - проекциям Р(х,у), по- 
лученным в результате измерений вне объекта [4-6]. 

СО ОО 

Р(р,Ѳ) = [ ] ст(х,т)-(5(/і-Т'С08(Ѳ)+Х-8Іп(Ѳ))с/хф', (1) 

—00 —00 

где 5(х,у) - фильтрующая дельта-функция Дирака. 


Решив ур. (1) относительно о(х,у), получаем 
формулу, позволяющую восстанавливать распреде- 
ления зарядов по проекциям: 


а(х,у) = 


2л 1 О 


2л о 
[ сІѲ ] Лр- 


ер(ѳ,р)/8 Р 


р-у-со$(0)+х-$т(0) 


( 2 ) 



Рис. 1. Получение проекции Р(р,Ѳ) при фиксированных зна- 
чениях р и Ѳдля распределения зарядов о(х,у) на по- 
верхности проводника 


Интеграл (2) существует в смысле главного зна- 
чения [4, 6]. Смысл переменных в формулах (1) и (2) 
ясен из рис. 1. В электростатическом поле потен- 
циал Щх,у) проводника в любой его точке одина- 
ков, заряд же распределен по поверхности неравно- 
мерно сфх,у). Поэтому основной задачей при опре- 
делении емкости является определение распределе- 
ния зарядов по поверхности проводника при задан- 
ном потенциале проводника. Чтобы получить связь 
между зарядом проводника и потенциалом на его 
поверхности необходимо осуществить обратное 
проецирование. Суть обратного проецирования в 
следующем. Проекция Р(р,Ѳ) двумерной функции 
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о(х,у) является одномерной функцией для каждого 
значения угла Ѳ. Если просуммировать исходный 
набор проекций Р(р,Ѳ) по углу Ѳ для §<Ѳ<2л, полу- 
чим суммарное изображение: 

*Дг 0 ) = — І СТ(Г) ,Л% г = {х,у},г 0 ={х 0 ,у 0 }. (3) 
4 п І г-г оІ 

Здесь 

2 п 

щт 0 ) = \р{р,ѳ)аѳ. ( 4 ) 

о 

Соотношение (3) в томографии известно как об- 
ратное проецирование ( Ъаск ргсуесгіоп) [4, 6], оно по- 
казывает, что распределение потенциала на поверхно- 
сти проводника есть низкочастотная версия распреде- 
ления зарядов на поверхности проводника [4, 5]. 

Авторами работы разработан томографичекий 
пакет программ, позволяющий восстанавливать 
распределение зарядов по его проекциям, а также 
рассчитывать емкости плоских пластин произволь- 
ной конфигурации. 

В качестве примера работы предлагаемого авто- 
рами алгоритма рассмотрим модельную задачу вос- 
становления распределения зарядов плоского 
электрода произвольной конфигурации. 

Сначала зададимся произвольным распределе- 
нием зарядов о{х,у) на плоской пластине радиусом 
К= 0,3 мм, порождаемым суперпозицией внешнего 
поля и поля, создаваемого самой пластиной, и опре- 
делим проекции этого распределения. Это есть ре- 


шение прямой задачи. Затем решим обратную зада- 
чу - по заданным проекциям восстановим распре- 
деления зарядов на поверхности электрода. При ре- 
шении обратной задачи необходимо учитывать осо- 
бенности решения обратных задач, т.к. такие задачи 
относятся к классу некорректных задач [6,7] и очень 
чувствительны к шумам в исходных данных. 

Пусть задано распределение зарядов на поверх- 
ности электродов: 


Ех(х -0.1, у-0.1, 30) 6 + 


<у(х,у) = \ 


+ Ех(х + 0Л,у-0Л,30) 6 - 
- Ех(х - Л ,у + ОЛ ,30) ь , г (х,у) <К 


0, г {х,у)>К 


где Ех(х,у,а)=ехр(-(х 2 +у 2 )сі:). Пространственный 
вид распределения зарядов о(х,у) приведён на рис. 2. 

Число дискретизаций области распределения 
зарядов по осям хи у выбираем равными N=52. Чи- 
сло дискретизаций области проекций по углу Ѳ и 
по параметру р выбираем равными соответственно 
ІѴ Ѳ = 1 8 , N=52. После использования дискретного 
аналога преобразование Радона - ур. ( 1), получаем 
распределение проекций изображения, рис. 3. 

Теперь решим обратную задачу. Используя вы- 
ражение (3), получаем распределение потенциала 
на поверхности пластины, вызванное внешним по- 
лем и полем, создаваемым самой пластиной. Решив 
интегральное ур. (3) путем сведения его к системе 
линейных алгебраических уравнений, находим 



Рис. 2. Пространственное распределение зарядов пластины о{х,у) и профили о(х,0) по оси х и о(0,у) по у соответственно 



Рис. 3. Пространственное распределение проекций Р(р, Ѳ) 
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Рис. 4. Результат реконструкции распределения зарядов: а) восстановленное распределение, 6) модельное распределение, 
в) пунктирные линии ~ модельные профили, сплошные ~ восстановленные 


о(х, у) -> ст = {от = а (г,),сг 2 = 

= а(г 2 ),ст 3 =а(г 3 ),...,ст л , = а(г іѵ )}: ( 5 ) 

ст = 4тг ■ А -1 • 1Т 0 , 

где А и =К( г ; ,г 0; ) - ядро интегрального ур. (3), 
і^еО...N - число точек дискретизации пластины. 

Выражение (5) является некорректным алгебра- 
ическим уравнением. Поэтому, осуществляя сла- 
бый спектральный сдвиг [4, 6-8], получаем заряд в 
дискретном множестве точек (рис. 4): 

сг = 4тг(А-А • Е) _1 ІІ 0 , (6) 

где Е - единичная матрица; Я=10 6 - параметр ре- 
гуляризации [2, 4-6]. 

При выбранных числах дискретизации ошибка 
восстановления распределения зарядов в каждой 
точке пластины не превышала 10 %. 

Приведём томографический метод расчёта ем- 
кости проводника в форме круглого кольца. Отно- 
шение внешнего радиуса кольца к внутреннему 
равно а=К/г. В осесимметричном случае преобра- 
зование Радона (1) преобразуется в преобразование 
Аббеля [4, 6, 8] 


Р(Р,Ѳ) = \^С*>°(х,у)х 


х8(р - у ■ С05(Ѳ) + х ■ сіхсіу — > 


->др) = 


2 | а(г)гйг 
\р\\/г 2 -р 2 


( 1 ) 


Интегрируя выражение (1) по углу, получим 
суммарное изображение - это есть распределение 
потенциала Щр)=2л-Р(р) на поверхности плоского 
электрода, рис. 5. 

Для вычисления заряда проводника будем ис- 
пользовать соотношение (2), которое в осесимме- 
тричном случае преобразуется в обратное преобра- 
зование Аббеля: 


а(х,у) 


[2 7Г СО 

у- 1 4Ѳ \ йр 

2л о -да 


дР(Ѳ,р)/8р 

р-у- СО8(0)+Х-8Іп(0) 


-> 


-»сг(р) = 


1 | Ц'(г)4г 

2л1 ІРІ \1>' 2 -р 2 


( 2 ) 


Теперь можно определить емкость проводника, 
используя известную формулу [1]: 




Рис. 5. Слева ~ распределение потенциала на поверхности проводника в форме круглого кольца, справа ~ его профиль 
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Рис. 6. Слева ~ фрагмент распределения заряда на поверхности проводника в форме круглого кольца, справа - его профиль 


С 


|| о(х,у)сіхсІу 

Е) 


V 


{х,ѵ}е/Х 


( 5 ) 


Для отношений величин внешнего и внутрен- 
него радиусов кольца сі=К/г, равного а~2, а 2 = 3, 
а 3 = 4, получены значения емкостей соответственно: 

С\ = 1,ШеК, С2=7,849аЯ, СЗ = 7,872еД 

Авторами были подсчитаны величины ёмкостей 
как канонических, так и неканонических форм. 
Приведём результаты расчетов ёмкостей различных 
пластин канонической формы на основе томогра- 
фического подхода и сравним их с расчетами [1]. 
Геометрические размеры пластин показаны на 
рис. 5-8. После определения распределения зарядов 
на поверхности электродов томографическим мето- 
дом, ёмкости пластин определялись по формуле (5). 

Результаты проведённых расчётов ёмкостей 
представлены на рис. 7-10. На рисунках приведе- 
ны верхние и нижние границы значений емкостей 
[1], а по центру проходят кривые, рассчитанные то- 
мографическим методом. 
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Рис. 7. Емкость пластины в форме ромба 



Рис. 8. Емкость прямоугольной пластины 



Рис. 9. Емкость пластины в форме кругового сегмента 
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Рис. 10. Емкость пластины в форме кругового сектора 

Результаты расчетов ёмкостей некоторых пла- 
стин неканонической формы на основе томогра- 
фического подхода приведены на рис. 11. Распре- 
деление потенциала на поверхностях проводников 
- величина постоянная, она повторяет форму 
электродов (рис. 11,6), тогда как заряды на поверх- 
ности электродов распределены неравномерно 
(рис. 11, в). 



Заключение 

Разработан томографический метод расчета ра- 
спределения зарядов на поверхности плоских про- 
водников. При известном распределении зарядов 
можно реконструировать электростатические поля 
и вычислять ёмкости проводников любой формы. 
Преимущество томографического подхода заклю- 
чается в возможности подсчитывать распределения 
зарядов, напряжённостей и емкостей проводников 
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г - С = 5.794-Я-е С = 6,215-К-е С=5.524-К-е 

Рис. 11. Результаты расчетов ёмкостей пластин неканонической формы: а) геометрические размеры, б) распределение потен- 
циалов и в) зарядов на поверхностях, г) соответствующие величины ёмкостей проводников 


неканонической формы. Современный уровень ра- 
звития компьютерной техники обеспечивает высо- 
кую точность расчетов. На основе предлагаемого 
подхода возможно решение обратной задачи мето- 
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